Topological properties of Sb(lll) surface: A density functional theory study 
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By using first-principles plane wave calculations, we systematically study the electronic proper- 
ties of the thin film of antimony in (111) orientation. Considering the spin-orbit interaction, for 
stoichiometric surface, the topological states keep robust for six bilayers, and can be recovered in 
the three bilayer film, which are guarantied by time-reversal symmetry and inverse symmetry. For 
reduced surface doped by Bi or Mn atom, localized 3-fold symmetric features can be identified. 
Moreover, the non-trivial topological states stand for non-magnetic substituted Bi atom, while can 
be eliminated by adsorbed or substituted magnetic Mn atom. 
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I. INTRODUCTION 

A new state of quantum matter-topological insulator 
(Tl)-has recently attracted great interest in condensed 
matter physics^. The realization of TI HgTc both from 
theoretical prediction and experimental observation 3 -^, 
open up opportunities for its potential application in 
semiconductor spintronics. It is believed that the quan- 
tum spin Hall effect as well as the time-reversal symmetry 
plays an important role in the new material, protecting 
the system from being disturbed by small perturbation 
caused by defects. This urges people to desire for more 
other promising materials, such as the Bi-based alloy 
Bii_ x Sb x , and layered compound Bi 2 X 3 (X=Se, Te) 5 - - — . 
Recently, the semimetallic antimony has been drawing 
lots of attention by its novel topological properties, and 
has become the proto-type system of Tl2r— . 

Distinct from Bi, which is topologically trivial, accord- 
ing to the parity of the band at the T point 13 , Sb has 
been identified to be a strong Tlii. Experimentally, it 
has been pointed out that the spin-split surface bands of 
Sb within its bulk band gap are connected to the con- 
duction band and valence band 9 -^. Moreover, the topo- 
logically nontrivial Sb thin films exhibit novel properties 
and provide a promising basis for spintronics applications 
such as device design and integration^ 2 .. 

Nevertheless, it is apparent that the detailed analysis 
of the electronic structures of Sb is lacking, and there 
are still many unanswered questions. Especially for the 
Sb(lll) surface, it is much desirable to identify the layer- 
dependence of the topological states of Sb thin film, and 
the modulating of the surface properties by impurities 
with or without magnetic momen t 15 ' 16 . Therefore, it is 
instructive to investigate the electronic structures of Sb 
to deeply understand its topological properties. This is 
not only of fundamental conceptual interest for under- 
standing the properties of TIs, but also paves the way 
for realizing promising application of the new state of 



matter. 

In this paper, we systematically study the properties of 
Sb(lll) surface by means of first-principles calculations. 
The bulk band structure of Sb is presented to reveal the 
domination of spin-orbit interaction (SOI) for topological 
properties of Sb. For the stoichiometric Sb(lll) surface, 
we investigate the layer-dependence of the surface states 
of Sb thin film. Moreover, we calculate the surface prop- 
erties of Sb when impurities arc introduced, including 
substituted nonmagnetic Bi and magnetic Mn, as well as 
the adsorbed Mn, and interesting results are obtained. 

The rest of the paper is organized as follows. In Sec. 
II the computational methods are briefly described. In 
Sec. Ill we present and discuss our results for the bulk 
and surface properties of Sb. Finally, in Sec. IV, we close 
our paper with conclusions of our main results. 



II. COMPUTATIONAL METHODS 

The calculations are performed within density func- 
tional theory using the Vienna ab-initio simulation pack- 
age (VASP)i£. The PBEi£ generalized gradient approx- 
imation and the projector-augmented wave potential 19 
are employed to describe the exchange-correlation energy 
and the electron-ion interaction, respectively. Here the 
Sb 5s and hp electrons are treated as valence electrons. 
The SOI that has been confirmed to play an important 
role in the electronic structure of TI, is included during 
the calculation. The cutoff energy for the plane wave ex- 
pansion is set to 400 eV. And a Fermi broadening 2 ^ of 
0.1 eV is chosen to smear the occupation of the bands 
around the Fermi level by a finite- T Fermi function and 
extrapolating to T = K. 

The Sb(lll) surface is modeled by a slab composing 
of several(l-6) bilayers (BL), and a vacuum region of 20 
A. Obviously, existence of the surface breaks the equiva- 
lence of the Sb atomic positions in a layer, hence through- 
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FIG. 1: (Color online) Atomic structure of bulk Sb and 
Sb(lll) surface, (a) Primitive cell of bulk Sb, (b) hexago- 
nal unit cell of bulk Sb, (c) slab model of Sb(lll) surface. 



out this paper, we employ the notation Sbl for the first 
atomic plane, and Sb2 for the second atomic plane from 
the surface, which are shown in Fig. 1(c). Integration 
over the Brillouin zone is done using the Monkhorst-Pack 
scheme— with 21x21x1 grid points for the p(lxl) sur- 
face; for the p(3x3) surface, 7x7x1 grid points are used. 
The structures of slab are fully optimized until the max- 
imum residual ionic force is below 0.02 eV/A. 



III. RESULTS AND DISCUSSION 




FIG. 2: Bulk band structure without SOI (a) and with SOI 
(b). The Fermi level is set to zero. 
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FIG. 3: Surface band structure with one to six bilayer thick- 
ness. The Fermi level is set to zero. 



The crystal structure of Sb is rhombohcdral with the 
space group Z)| d (i?3m), with two bismuth atoms in the 
trigonal unit cell (see Fig. Ha)). Equivalently, it can 
be represented in terms of a hexagonally arranged layer 
structure, which is shown in Fig. Hfb). The hexagonal 
unit cell has three sets of bilayers, where each bilayer con- 
sists of two Sb atoms. Structurally, bilayers in Sb form a 
stable unit with strong intrabilayer bonds, while the in- 
terbilayer bonding is much weaker. The calculated struc- 
ture parameters are a — 4.392 A, c = 11.432 A, which are 
close to experimental data a = 4.3007 A, c = 11.222 A22. 
Moreover, the optimized internal parameter is 0.233, con- 
sistent well with the experimental measurement 0.234. 

It has been proposed that Sb is a strong TI by calcu- 
lating its Z2 invariant v (=1) from the knowledge of the 
parity of the occupied Bloch wave function at the time- 
reversal invariant T point in the Brillouin zone 14 . This 
means that the SOI dominates the electronic nature of 
Sb, which can be seen from the band structure (see Fig. 
0). It is clear that the band gap is very small without 



SOI. While as the SOI is taken into consideration, the 
band gap is enlarged up to about 0.237 eV. Moreover, 
because of its topological electronic nature, we can ex- 
pect the existence of gapless surface states. Therefore, in 
the following we will focus our attention on the properties 
of Sb(lll) surface. 

It is well known that the topological features are of 
strongly thickness dependence of the TI thin film. For 
thinner film, the coupling between top and bottom sur- 
faces is strong enough to open up a whole insulating gap. 
While with increasing thickness, the inter-surface cou- 
pling becomes weaker and the topological features will 
be recovered. This is the same case for Sb(lll) surface. 
As illustrated in Fig. [31 we show the evolution of band 
structure of Sb films with the thickness from the sin- 
gle BL to six BL. In a single BL film, the Sb electronic 
states (mainly from 5p orbital) split into two parts form- 
ing a gap around the Fermi level. The gap is as large 
as about 0.9 eV, implying the strong coupling between 
two surfaces. In the case of two BL film, the splitting of 
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FIG. 4: (Color online) Surface charge density distributions at 
the height of 2.0 A above the reduced Sb(lll) surface: (a) 
Bi-substituted, (b) Mn-substituted and (c) Mn-adsorbed. 




5p states decrease with the decline of the coupling, lead- 
ing to a semimetallic electronic structure. Obviously, the 
topological features start to appear in the three BL case, 
where a non-trivial helical edge state (^o=l) below the 
Fermi level at T point can be identified. This helical 
state is consist of two surface states degenerated at T 
point but separated in energy elsewhere by SOI^.. Nev- 
ertheless, it is noticeable that for the case of four BL and 
five BL films, the gaps are opened up again at T point. 
This can be attributed to the inverse asymmetry of the 
films with four and five BL. While the topological states 
should be recovered when the number of BL of the film 
is multiples of three, which conserves the inverse sym- 
metry. As expected, it can be seen that the topological 
states are recovered for six BL. The double degenerate 
Sb(lll) surface states contain a single Dirac cone at the 
r point, which is robust and topologically protected by 
time-reversal symmetry. The Dirac point is about 160 
meV below the Fermi level, within the bulk band gap. 
Compared with the experimental value of 230 meV— , the 
difference may arise from the subsurface defects observed 
in the experiment. 

Various impurities (magnetic or not) may have differ- 
ent impact on the degeneracy and topological properties 
of Sb energy bands, which will be illustrated in the fol- 
lowing discussions based on the six BL film. It has been 
shown that the alloy Bii^Sb^ is a 3D Tli^22i2£. Here 
it would be interesting to investigate the surface state 
properties of reduced Sb(lll) surface with the existence 
of Bi as a nonmagnetic impurity. Moreover, as a compar- 
ison, a magnetic impurity (Mn atom) is also taken into 
consideration. During the calculations, the p(3x3) sur- 
face is adopted with the substituted impurities Bi or Mn 
atom taking the position of one subsurface Sb2 atom on 
each surface, while the adsorbed Mn atoms are symmet- 
rically located at the stable hep sites above two surfaces 
of Sb(lll). 

Figure 0] presents the surface charge density distribu- 
tions as simulation of STM topography at the height of 
2.0 A above the Bi-substituted, Mn-substituted and Mn- 
adsorbed surface, respectively. One can see clearly that 
the localized 3-fold symmetric features can be identified 



FIG. 5: (Color online) The band structure of clean (a) and 
Bi substituted (b) Sb(lll) surface. 



from all of three surfaces, especially for the Mn-doped 
ones. For the Bi-substituted surface, only tiny difference 
exists between the region above the Bi atom and others, 
and this can be attributed to that Bi and Sb belong to 
the same chemical group and possess the same number of 
valence electrons. As a result the Bi-doped Sb(lll) sur- 
face may maintain the topological features, which will be 
illustrated later by its band structure. For the Mn-doped 
surfaces, however, more evident different features can be 
observed. It is noticeable that the charge distribution 
is strongly depleted at the position just above the sub- 
stituted Mn atom, while the adsorbed Mn atom can be 
identified clearly by the charge accumulation. Keeping in 
mind that the magnetism of Mn atom, we will see that 
the magnetic impurity can have significant effect on the 
surface states of Sb(lll) surface. 

The band structure of Bi-doped Sb(lll) surface is pre- 
sented in Fig. EJb). For comparison, the band structure 
of stoichiometric p(3x3) Sb(lll) surface is also shown in 
Fig. EJa). It is clear that the topological surface state re- 
mains to be robust despite of the existence of Bi impurity, 
and the Dirac point of this reduced surface almost stays 
at the same position as that of the stoichiometric sur- 
face. This can be attributed to the following reasons, (i) 
The time-reversal symmetry stands against nonmagnetic 
impurities, (ii) The energy bands of Bi-doped system re- 
main degenerate thus keep gapless. Moreover, from the 
PDOS of Bi we can see that the doped Bi atom con- 
tributes little to the band structure near the Fermi level 
and the Dirac point, thus the topological features of Sb 
keep robust, which also consists with the fact that the 
alloy Bii^Sbz is a 3D TI inherited from Shii. 

Theoretically, if the impurity carries a magnetic mo- 
ment, the time-reversal symmetry is explicitly broken, 
and a local energy gap will be opened up at the Dirac 
point 1 ''. While up to now relevant experimental observa- 
tions are lack yet, therefore it is desirable to have a direct 
and intuitional sight into it. By calculations we find that 
the magnetic impurities can obviously lift the degener- 
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FIG. 6: The band structure of Mn-substituted Sb(lll) sur- 
face: (a) spin-up, (b) spin-down and (c) the PDOS of the Mn 
atomic d-bands. 



acy protected by the time-reversal symmetry. For clarity, 
here we plot the lifted band structures of Mn-substituted 
and adsorbed Sb(lll) surfaces dividedly in (a) and (b) 
of Fig. [5] and Fig. labeled by spin-up and spin-down, 
respectively. It can be seen that both of the magnetic im- 
purities eliminate the Dirac point by opening up a gap at 
T point. From the PDOS (See Fig. EJc) and Fig. [?_»), 
we find that unlike the doped Bi atom, the magnetic d 
orbital of Mn atom contributes much to the states near 
the Dirac point, hence clearly breaks the time- reversal 
symmetry. As expected^ a ferromagnetic ground state 
is formed on the TI surface by the introduced magnetic 
Mn impurity. 




FIG. 7: The band structure of Mn-adsorbed Sb(lll) sur- 
face: (a) spin-up, (b) spin-down and (c) the PDOS of the 
Mn atomic d-bands. 



IV. CONCLUSIONS 

In conclusion, we have systematically studied the topo- 
logical properties of Sb(lll) surface by density functional 
calculation. We found that the stoichiometric Sb(lll) 
surface possesses single Dirac point protected by the 
time-reversal symmetry and inverse symmetry. And the 
topological states are layer dependent and keep robust 
for six bilayers film. Moreover, we revealed that the non- 
trivial topological states stand for non-magnetic substi- 
tuted Bi, while the substituted or adsorbed magnetic Mn 
atom can obviously eliminate the Dirac point. We expect 
the present results for the topological features of Sb are 
greatly helpful for understanding and application of the 
topological insulator. 
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